Matter parity as the origin of scalar Dark Matter 



Mario Kadastik, 1 Kristjan Kannike, 1 and Martti Raidal 1 
1 National Institute of Chemical Physics and Biophysics, Ravala 10, Tallinn 10143, Estonia 

We extend the concept of matter parity Pm — (— l) 3 (- B-L > to non-supersymmetric theories and 
argue that Pm is the natural explanation to the existence of Dark Matter of the Universe. We 
show that the non-supersymmetric Dark Matter must be contained in scalar 16 representation (s) of 
SO(10), thus the unique low energy Dark Matter candidates are PM-odd complex scalar singlet(s) 
S and inert scalar doublet (s) H2- We have calculated the thermal relic Dark Matter abundance of 
the model and shown that its minimal form may be testable at LHC via the SM Higgs boson decays 
Hi -v DM DM. The PAMELA anomaly can be explained with the decays DM -> vlW induced 
via seesaw-like operator which is additionally suppressed by Planck scale. Because the SM fermions 
are odd under matter parity too, the DM sector is just our scalar relative. 



Introduction. While the existence of Dark Mat- 
ter (DM) of the Universe is now established without 
doubt pQ, its origin, nature and properties remain ob- 
scured. Any well motivated theory beyond the standard 
model (SM) must explain what constitutes the DM and 
why those DM particles are stable. In most popular mod- 
els beyond the SM, such as the minimal supersymmet- 
ric SM, additional discrete Z 2 symmetry is imposed by 
hand to ensure the stability of the lightest Z2-odd par- 
ticle. There is no known general physics principle for 
the origin of DM which could discriminate between the 
proposed DM models. 

In this Letter we propose that there actually might 
exist such a common physics principle for the theories 
of DM. It follows from the underlying unified symme- 
try group for all matter fields in grand unified theories 
(GUTs) and does not require supersymmetry. One can 
classify all matter fields in Nature under the discrete rem- 
nant of the GUT symmetry group which is nothing but 
the matter parity Pm- Thus the existence of DM might 
be a general property of Nature rather than an accidental 
outcome of some particular model. As a general result, 
there is no "dark world" decoupled from us, rather we 
are part of it as the SM fermions are also odd under the 
matter parity Pm- 

We argue that, assuming 50(10) [2] to be the GUT 
symmetry group, the discrete center Z n of U(l)x € 
50(10) remains unbroken. For the simplest case, n = 2, 
the GUT symmetry breaking chain 50(10) -> SU(5) x 
Pm implies that all the fermion and scalar fields of the 
GUT theory, including the SM particles plus the right- 
handed neutrinos Ni, carry well defined discrete quan- 
tum numbers which are uniquely determined by their 
original representation under 5O(10). We show that non- 
supersymmetric DM candidates can come only from 16 
scalar representations of 5O(10), and the unique low en- 
ergy DM fields are new SU(2) L x U(l)y Pj\/-odd scalar 
doublet (s) H 2 and singlet (s) 5 HE]. 

We formulate and study the minimal matter parity in- 
duced phenomcnological DM model which contains one 
inert doublet H 2 and one complex singlet 5. We show 



that the observed DM thermal freeze-out abundance can 
be achieved for wide range of model parameters. We also 
show that the PAMELA [5] and ATIC [7] anomalies in 
e + /(e~ + e + ) and e~ + e + cosmic ray fluxes can be ex- 
plained by DM decays via d = 6 [8] operators. In our case 
the Planck scale suppressed PM-violating seesaw-like op- 
erator is of the form m/(AxMp)LLH 1 H 2 , where m/Mp 
is Pjvf-violating heavy neutrino mixing. In this model the 
SM Higgs boson H x is the portal [9 to the DM. We show 
that for well motivated model parameter the DM abun- 
dance predicts the decay H\ — > DM DM, which allows 
to test the model at LHC [TO] . 

Matter parity as the origin of DM. The prediction 
of 5O(10) GUT is that the fermions of every generation 
form one 5O(10) multiplet 16i, i = 1,2,3. This is in a 
perfect agreement with experimental data as there exist 
15 SM fermions per generation plus right-handed N for 
the seesaw mechanism [llj . Assuming 5O(10) GUT, the 
first step in the group theoretic branching rule for the 
GUT symmetry breaking, 

5O(10) -> 5E/(5) x U(l) x -> SU(5) x Z 2 , (1) 

implies that every SU(5) matter multiplet [T5] and N 
carry an additional uniquely defined quantum number 
under the U(l)x symmetry. The U(l)x symmetry can 
be further broken to its discrete subgroup Z n by an order 
parameter carrying n charges of X |13[ 114] . The simplest 
case Z 2 , which allows for the seesaw mechanism induced 
by the heavy neutrinos N |llj . yields the new parity 
Px with the field transformation $ — > ±$. Therefore, at 
the electroweak scale after SU (5) symmetry breaking, the 
actual SM symmetry group becomes SU{2)l x U(1)y x 
Px- The discrete remnant of the GUT symmetry group, 
Px , implies the existence of stable DM. 

Under Pati-Salam charges B — L and the X-charge 
is decomposed as 

X = 'S(B-L)+AT 3R , (2) 

while the orthogonal combination, the SM hypercharge Y, 
is gauged in SU(5). Because X depends on 4T 3fl which 
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is always an even integer for T3R —1/2, 1, the Z 2 X- 
parity of a multiplet is determined by 3(B — L) mod 2. 
Therefore one can write 



Px=Pm = (-1? (B - L \ 



(3) 



and identify Px with the well known matter parity |15j . 
which is equivalent to i?-parity in supersymmetry. While 
U{l)xi X = 5(B — L)— 2Y, has been used to discuss and 
to forbid proton decay operators [TS], so far the parity 
([3]) has been associated only with SUSY phenomenology. 

Due to Eq. ([IJ a definite matter parity Pm is the gen- 
eral intrinsic property of every matter multiplet. The de- 
composition of 16 of 50(10) under Q is 16 = 1 16 (5) + 
5 16 (— 3) + 10 16 (1), where the U(l)x quantum numbers 
of the 5£/(5) fields are given in brackets. This implies 



= 16 



that under the matter parity all the fields 10 , 5 
l 16 are odd. At the same time, all other fields com- 
ing from small 50(10) representations, 10, 45, 54, 120 
and 126, are predicted to be even under Pm- Thus the 
SM fermions belonging to 16^ are all Pjvf-odd while the 
SM Higgs boson doublet is Pn/-even because it is em- 
bedded into 5 10 and/or 5 10 , and 10 = 5 10 (-2) + 5 10 (2). 
Although B — L is broken in nature by heavy neutrino 
Majorana masses, (— 1) 3 ( B ~ L ) is respected by interactions 
of all matter fields. 

As there is no DM candidate in the SM, we have to 
extend the particle content of the model by adding new 
5O(10) multiplcts. The choice is unique as only 16 con- 
tains Pjvf-odd particles. Adding a new fermion 16 is 
equivalent to adding a new generation, and this does 
not give DM. Thus we have only one possibility, the 
scalar(s) 16 of 5O(10). Because DM must be electrically 
neutral, 16 contains only two DM candidates. Under 
SU(2) L x U(l)y those are the complex singlet 5 = l 16 
and the inert doublet H2 € 5 16 . 

DM predictions of the minimal model. GUT 
symmetry groups are known to be very useful for classifi- 
cation of particle quantum numbers, and this is sufficient 
for predicting the DM candidates. Unfortunately GUTs 
fail, at least in their minimal form, to predict correctly 
coupling constants between matter fields. Therefore we 
cannot trust GUT model building for predicting details 
of DM phenomenology. Instead we study phenomenolog- 
ical low- energy Lagrangian for the SM Higgs Hi and the 
PM-odd scalars 5 and H%, 
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which respects Hi — > Hi and 5 — > —5, H 2 — > — H 2 . The 
doublet terms alone form the inert doublet model [5] . Fol- 
lowing Ref. [5], to ensure (5) = 0, we allow only the soft 
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Figure 1: Allowed 3<r regions for predominantly singlet DM 
in (ms, XsHi) plane for 6s = 5GeV, mn Q = 450 GeV. 



mass terms &g, ush a- n d the A5 term to break the inter- 
nal U(l) of the odd scalars. Thus the singlet terms in Q 
alone form the model A2 of [S] . The two models mix via 
^sh , Ush terms. Notice that mass-degenerate scalars are 
strongly constrained as DM candidates by direct searches 
for DM. The A5, 6| and ush terms in Eq. (j4| are crucial 
for lifting the mass degeneracies. 

We stress that our model of DM is based on the parti- 
cle quantum numbers and does not rely on numerology. 
However, the phenomenological studies of the model nec- 
essarily rise questions such as the gauge coupling unifica- 



tion. The one-loop /3-functions for gauge couplings g, g' 
and g 3 are given by f} g , = 7g' 3 , (3 g = -3g 3 and (3 g3 = 
—7g 3 . Based solely on the running due to those beta 
functions, we identify the unification scale 2 x 10 16 GeV 
by thesolution for g 2 = g 3 . The exact values of gauge 
couplings at M G are given by gi = y/5/3g' — 0.58, 
.92 = .93 = 0.53. The precision of unification of all three 
gauge couplings in our model is better than in the SM be- 
cause of the existence of an extra scalar doublet. We as- 
sume that an exact unification can be achieved due to the 
GUT thresholds corrections in full 5O(10) theory which 
we cannot estimate because the details of GUT symme- 
try breaking are not known [17 . In the minimal model 
with one extra doublet the required change of <?i at the 
GUT scale due to the threshold corrections is 10%. If, for 
example, there is one DM scalar multiplet for each gener- 
ation of fermions, the required threshold corrections are 
smaller, at the level of 4%. 

In the following we assume that DM is a thermal relic 
and calculate its abundance using MicrOMEGAs pack- 
age [H]. The DM interactions Q were calculated us- 
ing FeynRules package [19]. To present numerical exam- 
ples we fix the doublet parameters following Ref. [20 as 
m A ~ m H = 10 GeV, m H ± — mn = 50 GeV and treat 
TOjy and \i 2 as free parameters. For predominantly sin- 
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Figure 2: Allowed {mn , /12) parameter space for (j,sh = 
and different values of ms represented by the color code. 



glet DM we present in Fig.[T]the allowed 3er regions in the 
m s = Ms+^Sffi u2 /2— &| and \sh x plane for 65 = 5 GeV, 
m H a — 450 GeV and the values of fisn as indicated in 
the figure. For comparison we also plot the corresponding 
prediction of the real scalar model (light green band). For 
those parameters the observed DM abundance can be ob- 
tained for ms < mHo ■ Due to the mixing parameter /ish , 
a large region in the (m<j, X$Hi) plane becomes viable. 

To study DM dependence on doublet parameters we 
present in Fig.[2]the (mjj , ^2) parameter space for which 
the observed DM abundance can be obtained. Values of 
the singlet mass are presented by the colour code and we 
take ush — 0, 6s — 5 GeV. Without singlet 5, in the 
inert doublet model |20j . the allowed parameter space is 
the narrow region on the diagonal of Fig. [2] starting at 
m H ~ 670 GeV. In our model much larger parameter 
space becomes available. 

PAMELA, ATIC and FERMI data. PAMELA 
satellite has observed steep rise of e + /(e~ + e + ) cosmic 
ray flux with energy and no excess in p/p ratio [6]. ATIC 
experiment claims a peak in e~ + e + cosmic ray flux 
around 700 GeV [7], a claim that will be checked by 
FERMI satellite soon. To explain the cosmic e + excess 
with annihilating DM requires enhancement of the anni- 
hilation cross section by a factor 10 3-4 compared to what 
is predicted for a thermal relic. Non-observation of pho- 
tons associated with annihilation [5T] and the absence of 
hadronic annihilation modes |22 constrains this scenario 
very strongly. However, the PAMELA anomaly can also 
be explained with decaying thermal relic DM with life- 
time 10 26 s [25], 3-body decays in our case. 

In our scenario the global Z2 matter parity can be 
broken by Planck scale effects |T3] . If there exists, at 
Planck scale, a 50(10) fermion singlet N' , its mixing with 
the SU{5) PM-odd singlet neutrinos N via a mass term 
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Figure 3: Allowed 3<r regions in the singlet DM and SM Higgs 
boson mass plane for (is — and bs = 5 GeV. 

mNN' breaks Pm explicitly but softly. The exchange of 
N now induces also a seesaw-like [IT] operator 

TT~ T7~ LLHi H 2 IO- 30 GeV" 1 urW+H$, (5) 

where we have taken Ajv ~ 1, M/y ~ 10 14 GeV and m ~ 
v ~ 100 GeV. Such a small effective Yukawa coupling 
explains the long DM lifetime 10 26 s. 

LHC phenomenology. In our scenario the DM cou- 
ples to the SM only via the Higgs boson couplings Eq. Q . 
Therefore, discovering ~ 1 TeV DM particles at LHC is 
very challenging. However, if DM is relatively light the 
SM Higgs decays Hi — ► DM DM become kinematically 
allowed and the SM Higgs branching ratios are strongly 
affected. Such a scenario has been studied by LHC ex- 
periments 10] and can be used to discover light scalars. 

In our model such a scenario is realized for ^5 = 0, 
small bs <C v and heavy doublet. In this case the DM 
is predominantly split singlet and, in addition, the DM 
abundance relates the DM mass m| ~ Xsh^ 2 /2 — 6| to 
the SM Higgs boson mass , as seen in Fig. [3] For 
m Hl = 120 GeV, b s = 5 GeV we predict m_ s = 48 GeV 
with the Higgs branching ratios BR(Hi — > bb+cc+r?) — 
14.2%, BR(Hi -> DM DM) = 42.4% and BR{Hi -► 
S2 S2) — 42.4%. The second heaviest singlet ^2 with the 
mass m| 2 w Ash 1 w 2 /2 + b s decays via the SM Higgs 
exchange to S2 — > DM \i p, or S2 — > DM c c with almost 
equal branching ratios. Thus the SM Higgs boson decay 
modes are very strongly modified. This makes the Hi 
discovery more difficult at LHC but, on the other hand, 
allows the scenario to be tested via the Higgs portal 9 . 

Conclusions. We have extended the concept of Z 2 
matter parity, Pm = (—1) 3<,B ~ L \ to non-supersymmetric 
GUTs and argued that Pm gives the natural origin of 
DM of the Universe. Assuming that 5O(10) is the GUT 
symmetry group, the matter parity of all matter multi- 
plets is determined by their U(l)x charge under Eq. ([IJ. 
Consequently, the non-supersymmetric DM must be con- 
tained in the scalar representation 16 of 5O(10). This 
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implies that the theory of DM becomes completely pre- 
dictive and the only possible low energy DM candidates 
are the Pjvf-odd scalar singlet(s) S and doublct(s) Hi. 
We have calculated the DM abundances in the minimal 
DM model and shown that it has a chance to be tested at 
LHC via Higgs portal. Planck-suppressed Pm breaking 
effects may occur in the heavy neutrino sector leading to 
decays DM — > vlW which can explain the PAMELA and 
FERMI anomalies. 

Our main conclusion is that there is nothing unusual in 
the DM which is just scalar relative of the SM fermionic 
matter. Although B — L is broken in Nature by heavy 
neutrino Majorana masses, (— 1) 3 ( B ~ L ) is respected by in- 
teractions of all matter fields implying stable scalar DM. 
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